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Faraday patterns can be induced in Bose-Einstein condensates by a periodic modulation of the sys- 
tem nonlinearity. We show that these patterns are remarkably different in dipolar gases with a roton- 
maxon excitation spectrum. Whereas for non-dipolar gases the pattern size decreases monotonously 
with the driving frequency, patterns in dipolar gases present, even for shallow roton minima, a 
highly non trivial frequency dependence characterized by abrupt pattern size transitions, which are 
especially pronounced when the dipolar interaction is modulated. Faraday patterns constitute hence 
an optimal tool for revealing the onset of the roton minimum, a major key feature of dipolar gases. 



The formation of Faraday patterns in driven systems 
is a general phenomenon occurring in scenarios ranging 
from hydrodynamics and non-linear optics, to liquid crys- 
tals and chemical reactions l!| . Faraday patterns may be 
observed in Bose-Einstein condensates (BECs) by mod- 
ulating the nonlinearity arising from the interatomic in- 
teractions either by time-dependent Feshbach res- 
onances or by a time-dependent confinement. The 
latter method has been recently used for realizing these 
patterns in BECs [1]. Faraday patterns offer important 
insights about elementary excitations in BECs since the 
pattern size is determined by the Bogoliubov mode reso- 
nant with half of the driving frequency. For usual short- 
range interacting BECs the energy of elementary exci- 
tations grows monotonously as a function of their cor- 
responding momenta. As a result, the pattern size de- 
creases monotonously with the driving frequency. 

Recent experiments on atoms with large magnetic mo- 
ments [1, Q , polar molecules 0, Q , spinor Bose-Einstein 
condensates (BECs) ^ , and optical lattices are open- 
ing the rapidly-developing area of dipolar gases. In these 
gases, the dipole-dipole interactions (DDIs) play a signif- 
icant or even dominant role compared to the short-range 
interactions (SRIs). Dipolar BECs present a wealth of 



new physics ll|, [ij, [ij, [fj] due to the long-range and 



anisotropic character of the DDI. A major salient differ- 
ence between non-dipolar and dipolar BECs is provided 
by the dispersion of elementary excitations, which, due 
to the momentum dependence of the DDI, may show a 
roton-maxon character 14], similar to that encountered 
in Helium physics 15]. This key feature of dipolar gases 
may significantly influence, for deep roton minima, the 
BEC properties at finite temperatures [3] and even the 
BEC stability [13] • The roton minimum has not yet been 
observed experimentally, and it remains still an open 
question how to probe easily the onset of rotonization. 

In this Letter, we show that pattern formation is cru- 
cially modified in dipolar BECs with a roton-maxon spec- 
trum. Remarkably, contrary to many pattern forming 
systems, including non-dipolar BECs, the first unstable 
mode does not necessarily determine the emerging pat- 
tern, which may be dominated by harmonics of the driv- 



ing frequency with energies close to the roton minimum. 
As a result of that and of the multi-valued character of 
the roton-maxon spectrum the pattern size presents a 
highly non trivial dependence with the driving frequency 
characterized even for shallow roton minima by abrupt 
transitions in the pattern size. These transitions, which 
are especially pronounced when the DDI is modulated, 
may be employed to reveal easily the appearance of a 
roton- minimum in experiments on dipolar BECs. 

We consider a BEC of N particles with mass m and 
electric dipole d (the results are equally valid for mag- 
netic dipoles) oriented in the z-direction by a sufficiently 
large external field. The dipoles interact via a dipole- 
dipole potential: Vd{r) — — 3 cos^ {6)) / , where 9 
is the angle formed by the vector joining the interact- 
ing particles and the dipole orientation. We assume a 
strong harmonic confinement V{z) = muj'^z'^ /2 along the 
z-direction, whereas for simplicity of our discussion we 
consider no xy trapping. At sufficiently low tempera- 
tures the BEC wavefunction ^>{r) is given by the non- 
local non-linear Schrodinger equation (NLSE): 



(fr'U{f 



vl/(r), (1) 



where Hq = ~h'^\/^/2m + V{z) and U{r) = gS{r) + 
Vd{r — r'), with g = Airh^a/m, where a the s-wave scat- 
tering length (we consider in the following a < 0), and 
TO the particle mass. If the chemical potential (below 
introduced) ^2d f^z, the system can be considered 
"frozen" into the ground state (t)i^{z) of V and hence the 
BEC wave function factorizes as "^{r) — ^p{p)(j>o{z). Em- 
ploying this factorization, the convolution theorem, the 
Fourier transform of the dipole-potential and integrating 
over the z direction, we arrive at the 2D NLSE llSll: 
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where k is the a;?/-momentum, Iz = ^ h/mujz is the oscil- 
lator length, g2d = g/'J^h is the 2D short-range cou- 
pling constant, nik) is the Fourier transform of ]?/'(p)]^. 
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Figure 1: Dispersion of a 2D homogeneous BEC of ^^Ci with 
a — — 0.54nm (/3 — —0.375), a 3D density n^dj ^/^Iz ~ 
tujz — fJ,2d/2, and C, = h/y^Tmjl^i — 0.59fim. 
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Figure 2: Stability diagram (dark regions are stable) for the 
parameters of Fig. [T] as a function of a and k. (Left) hoj/^ = 
0.268 (regime B). (Right) hw/^i. = 0.134 (regime A). The most 
unstable modes are indicated by a solid line. 



and h2d(k) = 2 — erfc(fc), with erfc(fc) the com- 

plementary error function. The parameter /3 = (P/9 
characterizes the DDI strength compared to the SRL 

The homogeneous solution of ^ is tp{p,t) — 
\Jn2d exp [— z//2(iV^]i with nid the 2D density, and \i,2d — 
52d'T'2d(l + 87r/3/3) the chemical potential. The elemen- 
tary excitations of the homogeneous 2D BEC are plane 
waves with 2D wave number fc and dispersion [l^ 



,[kf = T(k) 



T{k) + 2g2dn2d 
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where T{k) = h^k^ /2m is the kinetic energy. If /3 = 
and since a < then e(fc ^ 0)^ < and phonon instabil- 
ity occurs, followed by the well-known collapse for attrac- 
tive short-range interacting BECs. This instability is pre- 
vented for sufficiently large DDI such that g+%T:(fi /?> > 0. 
At moderate d values, and due to the /c-dependence of 
the DDI (/i2d function), e(fc) may develop a roton-like 
minimum for intermediate k values (see Fig. [T]). The 
roton-maxon spectrum constitutes one of the most rele- 
vant novel features in dipolar gases. We show below that 
this roton minimum may be easily probed even for shal- 
low roton minima by modulating the system nonlinearity. 

We consider a modulation a{t) = ao[l -I- 2acos(2[ji)] 
about its mean cq, where a is the modulation amplitude 
and 2uj is the driving frequency. The homogeneous 2D so- 
lution is ■>pHip,t) = ^yn^exp[-i{t + ^ sin{2ujt))fj,2d/fi], 
with 7 = a/{l + SttP/S). The driving may induce a 
dynamical instability breaking the translational symme- 
try. This destabilization is best studied with an ansatz 
tp{p,t) = tpH{t)[l + w{t) cos(A:-r)], where is the com- 
plex perturbation amplitude. Inserting this ansatz into 
(12) we obtain a Mathieu equation for u =Ke{w): 



(Pu 1 . ^,^2 



2b{uj,k,a){hwfcos{2u;t)]u^0, (4) 



where b{uj,k,a) = 2a\g2d\n2dT {k) / {fajj^ , where g2d is 
calculated from the mean a^. Following Floquet Theo- 
rem, the solutions of dH) are of the form u[t) = c{t) exp at, 
where c{t) = c{t + 2-K/Ld) and cr(fc,a;,Q;) is the so-called 



Floquet exponent. If Re((T) > the homogeneous BEC 
is dynamically unstable against the formation of Faraday 
patterns, whose typical wavelength is dominated by the 
most unstable mode (that with the largest Re((T) > 0). 
For vanishing modulation the system becomes unstable 
at the parametric resonances e(fc) — nhco {n — 1, 2, ...). 

BECs with repulsive short-range interactions exhibit 
a spectrum e(fe)^ = T{k)[T{k) + 2g2dn2d], characterized 
by phonon-like excitations at low k and single-particle 
excitations at large k. For any given driving frequency 
the most unstable mode is always provided by the first 
resonance e(fc) = huj, and hence k — e~^(fta;) determines 
the typical inverse size of the Faraday pattern [23| ■ As a 
consequence a larger driving frequency leads to a pattern 
of smaller size, as shown in recent experiments In the 
following we show that the physics of Faraday patterns 
is remarkably much richer and involved in dipolar BECs. 

As mentioned above, for intermediate d values e{k) 
shows a roton minimum (with energy HcOr) and a maxon 
maximum [Hujm)- Hence, as a function of the modu- 
lation frequency 2uj we may distinguish three driving 
regimes: (A) a; < Wr, (B) LUr < uj < ujm and (C) 
Lu > LUm (see Fig. [1]). The latter regime is relatively 
uninteresting, since, as for non-dipolar BECs, the spec- 
trum is uni-valued and the most unstable mode is pro- 
vided by e(fc) = huj. The regime B on the contrary is 
multi-valued, and the condition e{k) = hiu is satisfied by 
a triplet ki < k2 < k^. These three resonant momenta 
lead to three instability tongues for growing modulation 
amplitude a (Fig. 2(a)). Hill's solution method 2^, 21 1 



provides that for the lowest resonance e{k) = fko, the 
Floquet exponent may be approximated for small b as 
(Ti = b{uj,k,a)/2 cx k"^ / {hwY . Hence the most unstable 
mode in regime B is always given by the largest momen- 
tum /ca, which dominates the Faraday pattern formation. 

For regimes B and C the Faraday pattern is, as for 
non-dipolar BECs, provided by e{k) — hw. The situation 
is remarkably different for regime A. The latter is better 
understood by considering the ratio (T2/cri between the 
Floquet exponents for the second (e(A:) = 2?iw) and the 
first (e(A:) = huj) resonance condition. This ratio may be 
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Figure 3: (Left) Ratio between the Floquet exponent cti for 
e — hu and CT2 for e = 2fia;. (Right) Most unstable fc as a 
function of uj. We use the same parameters as for Fig.[T] The 
dashed line indicates the roton frequency or momentum. 



Figure 4: Faraday patterns for hw/fi = 0.268 at t = 40.6 ms 
(left) and hu>/ jj, = 0.134 at i; = 93.2 ms (right), and the same 
parameters as in Fig. [1] 



obtained by using again Hill's solution method [13, Ef: pjg depicts the case of a frequency < luq < ^m, 
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where ^ = h/ ^2m^2d is the healing length. Not surpris- 
ingly, the first resonance dominates for a — ^ 0. However, 
contrary to the short-range interacting case, for a sur- 
passing a very small w-dependent critical a the situation 



changes completely. Fig. 3(a) depicts the ratio uijox as a 
function of lo for a small a = 0.04. Note that for a; > 
C2 < ci, and as expected, for regimes B and C the in- 
stability is dominated by the resonance e(fc) = hw. On 
the contrary for w < w^, cr2 > cri even for such a small 
value of a, and hence the lowest resonance is not any 
more the most unstable one. This surprising result is a 
direct consequence of the non-monotonous character of 
the roton-maxon dispersion law. 

Our numerical Floquet analysis shows indeed (see 
Fig. 3(b) I that for a > acr (for the parameters of Figs. [2] 
acr — 0.027) the most unstable mode for all driving fre- 
quencies within the regime A is given by the largest mo- 
menta k compatible with the first harmonic e(fc) = nhuj 
lying in the regime B (or, if none, the first lying in regime 
C). This has important consequences for the wavenum- 
ber selection as a function of the driving w, which, as 



shown in Fig. 3(b) is remarkably different than that for 



the non-dipolar case. The pattern size does not decrease 
monotonously with growing uj, but on the contrary oscil- 
lates in the vicinity of the roton momentum, presenting 
abrupt changes of the pattern size at specific driving fre- 
quencies. These oscillations are the result of the subse- 
quent destabilization of higher harmonics in regime B. 

This distorted wave number selection is directly mir- 
rored into the spatial form of the corresponding Faraday 
patterns. We have studied the dynamical instability in- 
duced by the modulation and the corresponding Faraday 
patterns by simulating Eq. ([2]) numerically with periodic 
boundary conditions and an overimposed random noise 
provided by a tiny random local phase (< 10 ^"^tt) on 
the homogeneous solution. Our direct numerical calcula- 
tions is in excellent agreement with our Floquet analysis. 



where as expected the Faraday pattern is indeed given by 
the largest resonant momentum. In Fig. |4(b)] we consider 
Lu = wo/2 which is within the regime A. Strikingly, due 
to the discussed selection of higher harmonics, the Fara- 
day pattern is basically the same as for a double driving 
frequency ujq. This quasi- insensitivity becomes quantita- 
tively evident after Fourier transforming the pattern. 

In the previous discussion we considered the modula- 
tion of the s-wave scattering length a{t). A dipolar BEG 
offers, however, an additional novel way of modifying the 
system nonlinearity by a time-dependent DDI. This may 
be achieved by modulating slightly the intensity of the 
polarizing field (e.g. the electric field orienting a polar 
molecule) or by introducing a slight precession of the di- 
rection of the external field (e.g. by additional transver- 
sal magnetic fields in the case of atomic dipoles). In the 
following we show that the Faraday patterns obtained by 
means of a modulated DDI differ very significantly from 
those obtained by modulating a{t). 

We consider a temporal modulation of the DDI (P — 
g/3{t), with f3{t) — f3[l+2a cos{2Lut)] about its mean value 
(3. Following a similar procedure as that discussed above 
for the case of modulated a{t), we obtain the Mathieu 
equation for the real part of the perturbation amplitude. 
This equation is of the same form as Eq. ([4]) but with 



b{u), fc, a) 



8TTa\Pg2d\n2d 



T{k)h2d [ ^ 



(6) 



The modified fc-dependence of h(LO, fc, a) has crucial con- 
sequences for the formation of Faraday patterns. Simi- 
larly as above we may obtain from Hill's solution method 
the Floquet exponent for the first resonance e(fc) = 
CTi (X k'^h2d{klz/V2). This leads to a remarkably dif- 
ferent selection rule for u; values within the regime B. 
Contrary to the case of modulated a{t), it is the interme- 
diate momentum fc2 and not the largest one fca the most 
unstable within regime B. This leads to a remarkably 
abrupt change in the Faraday pattern size in the vicinity 
of ujm [111 • In addition, and similar to the case of mod- 
ulated a{t), driving with uj < may be dominated by 
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Figure 5: Most unstable fc as a function of lj for modulated 
l3{t), a = 0.12, and the same parameters as for Fig. [1] We 
indicate the roton and maxon frequencies and momenta. 




Figure 6: Faraday patterns for huj/ = 0.32 at t = 116.7 ms 
(left) and huj/ fj, = 0.34 aX t = 405.1 ms (right), and the same 
parameters as in Fig. [5] 



higher harmonics. Fig. [5] shows the most unstable fc as a 
function of u) for a typical case of modulated (3{t). Note 
not only the above mentioned abrupt jump in the vicin- 
ity of but also at other lo values within the regime A. 
As for the case of modulated a{t) these jumps represent 
abrupt transitions in the Faraday pattern size, which are 
certainly much more marked than for the modulated a 
case. Figs. |6] show the abrupt change in the patterns for 
two driving frequencies right below and above the tran- 
sition close to Ldm [201 . 

In our calculations we have assumed for simplicity no 
trapping on the xy plane. An harmonic xy-confinement 
(with frequency ljj,) leads to a finite momentum cut-off 
kc = \/ mu!±/h. In a good approximation all features in 
the excitation spectrum with momenta k ^ kc are not af- 
fected by the inhomogeneous trapping. For typical roton 
momenta fc^ ~ 0.5 and C ~ O.Gfim in our figures, k^ kc 
demands for ^^Cr a transversal frequency uj±_ < 130Hz, 
which can be considered a typical experimental condition. 
Finally, we stress that Faraday patterns are a transient 



phenomenon, and that for the case discussed here (a < 0) 
pattern formation is followed by collapse (and consequent 
violation of the two-dimensional condition) . 

In summary, the physics of Faraday patterns is largely 
modified in dipolar BECs in the presence of even shal- 
low roton minima. Whereas in non-dipolar BECs the 
Faraday pattern size decreases monotonously with the 
driving frequency 2lj^ in dipolar BECs the patterns show 
a w-dependence characterized by abrupt changes in the 
pattern size, which are especially remarkable when the 
dipole itself is modulated. Faraday patterns constitute 
hence an excellent tool to probe the onset of rotonization 
in on-going experiments with dipolar condensates. 

This work was supported by the DFG (SFB407, 
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